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ABSTRACT 

We present a deep, wide-field optical study of the M81 group dwarf galaxy Holmberg II 
(HoII) based on Subaru/Suprime-Cam imaging. Individual stars are resolved down 
to I ~ 25.2, i.e. about 1.5 mag below the tip of the red giant branch (RGB). We 
use resolved star counts in the outskirts of the galaxy to measure the radial surface 
brightness profile down to \iy ~ 32 mag arcsec -2 , from which we determine a projected 
exponential scalelength of 0.70' ± 0.01' (i.e. 0.69 ± 0.01 kpc). The composite profile, 
ranging from the cored centre out to R=7', is best fit by an EFF profile which gives 
a half-light radius of 1.41' ± 0.04' (i.e. 1.39 ± 0.04 kpc), and an absolute magnitude 
My = —16.3. The low surface-brightness stellar component of HoII is regular and 
symmetric and has an extent much smaller than the vast H I cloud in which it is 
embedded. We compare the spatial distribution of the young, intermediate age, and 
old stellar populations, and find that the old RGB stars are significantly more centrally 
concentrated than the young stellar populations, contrary to what is observed in most 
dwarf galaxies of the Local Universe. We discuss these properties in the context of 
the comet-like distribution of H I gas around HoII, and argue for the presence of a hot 
intragroup medium in the vicinity of HoII to explain the contrasting morphologies of 
the gas and stars. 

Key words: galaxies: individual: Holmbergll - galaxies: dwarf - galaxies: irregular 
- galaxies: stellar content - galaxies: groups: individual: M81 group - intergalactic 
medium. 



1 INTRODUCTION 

In spite of their intrinsic faintness and minimal contribu- 
tion to the total light, the stellar outskirts of galaxies hold 
crucial information about the processes of galaxy formation 
and evolution. Interactions, mergers and accretions all leave 
their imprint on the outer stellar populations in the form of 
substructures, streams and diffuse haloes. In addition, the 
long evolutionary timescales and high sensitivity to external 
influences means that coherent substructures are easier to 
detect and longer-lived than in the denser inner regions. 
Considerable effort has been devoted to studying the 



* Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan, 
f E-mail: ejb@roe.ac.uk 



outer regions of massive galaxies (see, e.g. lBarker et al 120121 . 
and references therein) which, according to hierarchical 
models of galaxy formation, have acquired a significant frac- 
tion of their mass through mergers and accretion episodes. 
The stellar peripheries of low mass dwarf galaxies have been 
much less studied. Most work to dat e has focused on Lo- 
cal Group (LG) galaxies (e.g. Forn ax: IColeman et alj|2005l: 
NGC 6822: Ide Blok fc WalteiboOrj ; Sculptor: IWestfall et al l 
2006), although a few dwarfs in nearby groups have also 
been targeted fe.g. lRvs et alj201lf ). The strong morphology- 
density relation exhibited by dwarf galaxies suggests exter- 
nal mechanisms play a major role in shaping their evolution 
fe.g. lWeisz et alj|201ll 1; surveying their stellar outskirts may 
therefore yield clues on the dominant processes involved. 

The global picture emerging from previous work is that 
most dwarf galaxies have a smooth, and generally old and 
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Figure 1. Mosaic V image of HoII, showing the whole field-of- 
view o f our Subaru data (~35.8' X 29.3'). The Hi contours from 
lBC02t ranging from Njjj = 0.1 to 19 X 10 20 atoms cm 2 , are 
overplotted. Note the comet-like morphology of the outermost 
contour. 



metal-poor 'halo' surrounding a more concentrated, younger 
and/or more metal-rich population. While the stellar haloes 
of large galaxies are believed to have formed from the accre- 
tion of smaller galaxies at high redshift, it is unclear whether 
this process has been significant in dwarfs. For example, 
pure accretion cannot explain the main properties of dwarf 
galaxies' haloes, i.e., smooth distribution, and the existence 
of age and/or metallicity gradients. Instead, these obser- 
vations suggest either a 'shrinking' scenario, in which the 
region of active star formation contracts as gas supply di- 
minishes i|Hidalgo et al.ll2009l ; Izhang et al.ll2012l ), or radial 
migration of stars formed clo se to the centre towards the 
outskirts (jStinson et alJl2009l ). 

Here we analyse the stellar outskirts of the dwarf 
galaxy Holmbergll (hereafter abbreviated as HoII). HoII 
is a dwarf i rregular ga l axy in the M81 group that was dis- 
covered by iHolmbera l|l950l ) while surveying the galaxies 
in this group. Due to its location on the near-side of the 
group ((m— M)o=27.65, i.e. 3.4 Mpc), as well as its prox- 
imity to the Sc spiral galaxy NGC2403 and similar ra- 
dial velocity, it is usually associated with the NGC2403 
subgroup along with thr ee other dwarf irregular galaxies 
l)Karachentsev et alJl2003 ). HoII is very similar to the Small 
Magellanic Cloud (SMC) in terms of absolute magnitude, 
Hi and total mass: M s ~ -16.7, M H i ~ 6 x 10 s 
and Mtot = 2.1 x 1O 9 M ijWalter et all 120071 : lOh et all 
l201ll ) compared to M B ~ -16.1, M m ~4x 10 8 Mq and 
Mtot = 2.4 x 10 9 Mf7i for the SMC jStanimirovic et al]|l999t 
IStanimirovic. Stavelev-Smith. fc Jonesll20o3 )! 

While HoII has been observed at all wavelengths, to 
date the only deep resolved s tellar populations study comes 
from HST/ACS observations (Weisz et al. 2009) which cover 
a relatively small fraction of the galaxy. The need for 
deep wider data is especially motivated by the strik- 
ing morphology of its H i cloud. From deep VLA data, 
iBureau fc Carignanl l|2002l . hereafter BC02) found that the 
distribution of neutral hydrogen has a cometary appearance 
- compressed on one side with a faint extended component 



on the opposite side - wi th the tail pointing away from the 
centre of the M81 group. |BC02| argued that HoII is moving 
toward the M81 group and that ram pressure from a hot 
intragroup medium (IGM) is responsible for the Hi mor- 
phology, although they could not rule out the alternative 
interpretation of a gravitational interaction between HoII 
and one of its fainter neighbours. 

Possible signatures of ram pressure stripping have been 
obse rved in a number of galaxies of the local universe 
(e.g. ICondonlll983l: iGavazzi et al.lll995l ; iRvder et alJ I l997| : 
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However, most of these galaxies 
are too distant to be resolved into individual stars and their 
analyses have been confined to fairly high surface brightness 
inner regions. Knowledge of the distribution of resolved stel- 
lar populations at very large radii is fundamental because 
the stars do not respond to ram pressure. On the other hand, 
tidal forces affect gas and stars equally so both components 
should exhibit similar asymmetries. Comparing the large- 
scale distribution of the stars with that of the gas in HoII 
has the potential to reveal whether the H I morphology was 
caused by ram pressure, tidal forces, or a combination of 
both. 

In this paper we present a deep, wide-field study of HoII 
based on Subaru/SuprimeCam data, and analyse the prop- 
erties of its stellar populations in the context of the comet- 
like shape of the H I cloud in which it is embedded. In Sec- 
tion [21 we describe the observations and data reduction, and 
present the resulting CMDs in Section [3] The spatial dis- 
tribution of the various stellar populations are decribed in 
Section [3] In Section we present the radial profile and 
constrain the spatial extent of HoII. We discuss the impli- 
cations of our results regarding the peculiar stellar and H I 
distributions in Section [S] and summarize the main results 
in Section [7] 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations and Image Processing 

The observations of HoII (=UGC 4305 =DDO 50 =Arp 268) 
were o btained in service mod e with the Suprime-Cam instru- 
ment (|Mivazaki et all l2002f ) on the 8-m Subaru telescope 
(P.I.: M. Barker). The data were acquired on the night of 
December 18, 2009. Suprime-Cam is a mosaic camera made 
of 10 chips disposed in two rows of five chips - thus produc- 
ing 10 images per exposure - leading to a large field-of-view 
(FOV; 34' x27'). A single pointing was therefore sufficient to 
cover t he relatively small g alaxy (R25 = 4.1'; from HYPER- 
LEDA: iPaturel et al.ll2003l ). in cluding the large Hi cloud in 
which it is embedded (R ~ 16'; |BC02| ). This is illustrated in 
Figure [TJ showing the V-ban d mosa ic from our Subaru data 
where the H I contours from BC02 have been overlaid. 

Ten exposures were obtained in each band, with indi- 
vidual exposure times of 600s in Johnson V and 240s in 
Cousin /. A small dithering pattern was used to cover the 
gaps between the chips as well as limit the effects of bad 
pixels and other camera defects, resulting in a mosaic image 
of ~35.8' x29.3'. The observations were carried out in pho- 
tometric conditions, with seeing in the range 0.65 — 0.95". 
However, due to HoII's location at high declination, it can 
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only be observed at relatively large airmasses from Hawaii 
(>1.6), so the median seeing measured in the V and I images 
is 0.96 and 0.76", respectively. 

The i mage processing procedure s were similar to those 
followed bv lBarker et alj (|2009| . |2012| ). After converting each 
exposure to a single multi-extension FITS file, all images and 
calibration frames were run through a variant of the data re- 
duction pipeline developed for processing Wide Field Cam- 
era (WFC) data from the Isaac Newton T elescope (INT) 
- for father details see llrwinl (| 19851 . Il997l ); llrwin fc Lewis] 
|200ll ); llrwin et all (|2004T ). 

First stage image processing steps included bias and 
overscan-correction, together with trimming to the reliable 
active detector area. Master flats were created by stacking 
a dithered set of 10 V and 12 I twilight sky exposures. The 
flat-fielding step also corrects for internal gain variations 
between the detectors. After flat-fielding the dark sky I- 
band images were examined for signs of fringing, but as for 
other similar Subaru data, this was found to be negligible. 



Prior to stacking, detector-level catalogues were gener- 
ated for each individual processed science image to refine 
the astrometric calibration and also to assess the data qual- 
ity. For astrometric calibration, a Z enithal polynomial pro- 
jection (|Greisen k, Calabretta|[2002l ) was used to define the 
World Coordinate System (WCS). A fifth-order polynomial 
includes all the significant telescope radial field distortions 
leaving just a six-parameter linear model per detector to 
completely define the required astrometric transformations. 
The T wo Micron All Sky Survey (2MASS) point-source cat- 
alogue (|Cutri et al.l 12003 ) was used for the astrometric ref- 
erence system. 

During the stacking process the individual Subaru cat- 
alogues were used in addition to the standard WCS solu- 
tion, to further refine to the sub-pixel level the alignment of 
the component images. The common background regions in 
the overlap area from each image in the stack were used to 
compensate for sky variations during the exposure sequence 
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Figure 3. Completeness of the Subaru data, relative to the 
HST/ACS data, as a function of galactocentric radius for stars 
brighter than 7rj = 24.5. The dashed lines is a fit to the points 
between R=1.5 and 4'. 

and the final stack included seeing weighting, confidence (i.e. 
variance) map weighting and clipping of cosmic rays. 

As a final image processing step, catalogues were de- 
rived from the deep stacks for each detector and their WCS 
astrometry was updated prior to forming the mosaic over all 
detectors. Any residual small offsets in the underlying sky 
level between each detector stack were removed iteratively 
by visual inspection of a 4x4 blocked mosaic of the whole 
field. 

Full resolution mosaics were used to provide an initial 
set of full-fiel d catalogues base d on standard aperture pho- 
tometry (e.g. llrwin e t al. 2004j). Since no standard star field 
was observed, we used INT WFC V- and i-band observations 
centred on HoII and taken in photometric conditions during 
April 2009 to derive the photometric calibration for the mo- 
saics. The INT WFC V,i photom etry was first converted to 
the Johnson-Cousins system (e.g. iMcConnachie et al.ll2003T ) 
and then directly used to calibrate the Subaru data with an 
estimated systematic error of ±2 per cent. 

Figure [2] shows a colour composite image obtained from 
the V and I mosaics, which highlights the irregular morphol- 
ogy and strong contribution from the young, blue supergiant 
stars. 



2.2 Profile-fitting Photometry 

Given the very high stellar density in HoII, we subsequently 
decided to perform profile-fitting photometry which is bet- 
ter suited to obtain accurate measurements in crowded 
fields. The PSF photometry was carried out with the 
DAOP HOT / ALLSTAR/ ALLFRAME suite of programs dStetsonl 
1 1994 ) as follows. We performed a first source detection at 
the 3-cr level on the individual images, which was used as 
input for aperture photometry. From these catalogues, 300 
bright, non-saturated stars per image were initially selected 
as potential PSF stars. An automatic rejection based on the 
shape parameters was used to clean the lists, followed by a 
visual inspection of all the stars to remove the remaining 
unreliable stars. We ended up with clean lists containing at 



least 150 good PSF stars per image. Modelling of the em- 
pirical PSF with a radius of 12 pixels was done iteratively 
with DAOPHOT: the clean lists were used to remove all the 
stars from the images except PSF stars, so that accurate 
PSFs could be created from non-crowded stars. At each it- 
eration, the PSF was modelled more accurately and thus 
the neighbouring stars removed better. Every 3 to 5 iter- 
ations, the degree of PSF variability across the image was 
also increased, from constant to linear, then quadratically 
variable. 

The following step consisted of profile-fitting photome- 
try on the individual images using ALLSTAR with the empir- 
ical PSFs previously created. The resulting catalogues were 
matched on a chip-by-chip basis, keeping only the objects 
for which the PSF fitting converged in at least three images 
per band to limit the number of false detections, to create 
one clean stellar catalogue per chip. These catalogues were 
further cleaned by rejecting extended objects based on the 
sharpness parameter, then merged to produce a master cata- 
logue for the whole Suprime-Cam FOV. This catalogue was 
then used as the input star list for ALLFRAME, which was 
run on all the individual images at once. The output of all- 
frame consists of a catalogue of PSF photometry for each 
individual image. A robust mean magnitude was obtained 
for each star by matching these catalogues. 

The final photometry was calibrated to the Johnson- 
Cousin standard system by matching ~1750 bright stars in 
common with the calibrated aperture photometry described 
in the previous Section. The uncertainty on the offset be- 
tween the aperture and PSF photometry is smaller than 
0.001 in both bands and therefore represents a negligible 
contribution to the to tal magnitude uncertainties. 

According to the [Schlegel, Finkbeiner, fc Davis! (|l998l ) 
reddening maps, the area covered by our observations suffers 
from minor foreground differential reddening, with E(B— V) 
ranging from about 0.023 to 0.033. To obtain accurate pho- 
tometry over the whole FOV, each individual star was cor- 
rected for reddening based on its location. We note that the 
corresponding dust mask map indicates that a small area 
centred on HoII (out to r~6'; i.e. ~10 per cent of the Sub- 
aru FOV), as an extragalactic source, was removed from the 
dust map and replaced with the median value of the sur- 
rounding pixels. This means that the star-by-star correction 
does not include extinction internal to HoII. 

Finally, the astrometric calibration for the whole stel- 
lar catalogue was obtained with the IRAF tasks ccxymatch, 
ccmap, and cctran using ~800 stars in common with V ersion 
2.3.2 of the Guide Star Catalog II l|Lasker et alj|2008h . The 
accuracy of the resulting astrometry is about 0.3". 

2.3 Completeness 

To estimate the completeness of our data, we retrieved the 
HST/ACS photometry of HoII from the ACS Nearby Galax y 
Survey Treasury program (ANGST; Ibalcanton et alj|2009r i. 
While it only covers a tiny fraction of our FOV (~2 per 
cent), it is located on the highest density area which is the 
most affected by incompleteness. It is ~3 mag deeper than 
our Subaru photometry and much less affected by stellar 
crowding thanks to the higher spatial resolution of the in- 
strument. Therefore, in the following we assume that it is 
100 per cent complete in the range of magnitudes covered 
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by our Subaru data, and representative of the intrinsic pho- 
tometric properties of the stars in HoII. 

We used daomatch and daomaster (|Stetsonll 19931 ) to 
match the HST and Subaru photometric catalogues. These 
programs are based on the robust 'matching triangles' tech- 
nique, rather than assuming a given matching radius. This 
ensures a proper matching of the stars even in very crowded 
regions, regardless of the possible translations, rotations, 
scale changes, or flips of the coordinate systems. 

Using all the stars from the area in common between 
the two photometric catalogues, we find the 50 per cent 
completeness limits of the Subaru data at Vb=23.93 and 
Jo=23.00. However, the centre of the galaxy is more affected 
by crowding than the outer regions, so completeness varies 
significantly with radius. To correct the radial density profile 
shown in Section [S] which uses all the stars with Iq < 24.5, 
we therefore estimate the completeness down to this magni- 
tude as a function of radius. The result is shown in Figure [3] 
For each annulus, the errorbar was calculated by adding in 
quadrature the Poisson uncertainties on the number of stars 
in the HST and Subaru subsamples. We find that complete- 
ness varies from 10 per cent in the inner 0.5' to 90 per cent at 
R~4'. The dashed line is a fit to the points between R=1.5 
and 4', and suggests that we can assume ~100 per cent com- 
pleteness beyond 4'. 



3 COLOUR-MAGNITUDE DIAGRAMS 

The resulting colour-magnitude diagrams (CMDs) for the 
field and HoII stars are shown in Figure [4] We used the 
0.1 and 2xl0 20 atoms cm H I contours to separate the 
two populations (see Section [2} : stars outside the former are 
assumed to be field stars while stars inside the latter are 
considered HoII. The CMDs were cleaned of non-stellar ob- 
jects using the photometric parameters given by ALLFRAME, 
namely oy,i *S 0.2 and | SHARP | ^ 1. To enhance the fea- 
tures at faint magnitudes, the CMD of the top left panel 
was further cleaned using tighter constraints on the sharp- 
ness (|SHARP| ^ 0.5) but is shown only for illustrative pur- 
pose. It contains ~16,600 stars, instead of the ~27,100 stars 
in the CMD of the top right panel. 

At J brighter than ~22 the field CMD, shown in the 
bottom right panel of Figure[4] harbors two prominent verti- 
cal sequences: the main-sequence turn-off stars of the Milky 
Way (MW) halo at (V - J)o ~ 0.7, and dwarf stars in the 
MW disc at (V — I)o ~ 2.5. At fainter magnitudes, these se- 
quences are overwhelmed by the contribution of unresolved 
background galaxies. We note, however, that the HoII CMDs 
in the top panels correspond to an area about 10 times 
smaller than that covered by the field stars, so the effect 
of the contamination by foreground stars and background 
galaxies is not as significant. To illustrate this, the bottom 
left panel shows the CMD of a field with the same area as 
HoII. The scaled field CMD contains ~7600 objects. 

To help identify the features in the HoII CMDs, in the 
top right panel of Figure [3] we overplot isochrones from the 
Padua stellar evolution library l|Marigo et al.1120081 ). shifted 
to the distance of HoII. The young isochrones have Z=0.002 
and ages of 10, 20, 50, 100, and 160 Myr from top to bot- 
tom, while the old ones shown in red are 12.5 Gyr old with 
Z=0.0001, 0.0003, 0.001, and 0.002 from left to right. 




Figure 4. Extinction-corrected color-magnitude diagrams of 
HoII (top panels) and field stars (bottom panels). In 
the bottom right panel, the contour levels correspond to 
[6,12,18,24,30] xlO 3 stars mag , and the error bars show the 
mean photometric errors as a function of magnitude for both 
field and HoII stars. The boxes used to select the various stellar 
populations are shown in the right panels: MS+BSG (blue), RSG 
(yellow), AGB (red ), and RGB (orange ). Selected isochrones from 
the Padua library llMarigo et al.|[2008l) have been overplotted in 
the top right panel (see text for details). 



The most prominent features of the HoII CMD are the 
blue and red supergiant branches (BSC and RSG, respec- 
tively) of core helium-burning stars at (V — I)o ~ and 
1.2, respectively. Their presence indicates vigorous star for- 
mation in the past ~160 Myr. The BSG branch curves to- 
ward the red below Io ~ 24, and joins the RSG close to 
the red clump just below our limiting magnitude. Redward 
of the RSG and at Io > 23.5, a higher density of stars 
corresponding to the red giant branch (RGB) is visible. 
The expected magnitude of th e tip of the RGB (TRGB), 
from iKarachentsev et al.l (|2002l ) , is shown as a thin line at 
Io = 23.60. Some of the stars brighter than the TRGB and 
redder than the RSG are genuine asymptotic giant branch 
(AGB) stars (see Section 0J, although they are difficult to 
separate from the contamination of foreground stars and un- 
resolved galaxies. Finally, a few dozen main sequence (MS) 
stars forming a diffuse band blueward of the BSG branch 
are visible at about (V — J)o = —0.2 and Io > 21.5, in- 
dicating star formation within the last 10 Myr. We note 
that while some of these stars lie blueward of the youngest 
isochrones, their concentration close to the centre of HoII 
confirms that they are bona fide members of HoII. Their 
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Figure 5. Spatial distribution of candidate MS+BSG, RSG, AGB and RGB stars in the whole Subaru FOV (left Vand the corresponding 
contours (right). Stellar contours levels are arbitrary. The two Hi contours (0.1 and 2xl0 20 atoms cm -2 , from BC02) used to separate 
HoII and field stars are shown as thick gray lines to use as guides. 



color spread, larger than expected from photometric errors 
only, is a consequence of the crowding in this region. 

In the following Section, we analyse the spatial distribu- 
tions of the different stellar populations individually. They 
are selected on the CMDs using the boxes shown in the right 
panels of Figure [4] The blue, yellow, red, and orange boxes 
enclose the stars belonging to the MS+BSG (<160 Myr), 
RSG (10-160 Myr), AGB (a few hundred Myr to a few Gyr), 
and RGB (>1.5 Gyr), respectively. For the MS+BSG selec- 



tion, we chose to use the whole range of luminosity since 
contamination by foreground and background sources is vir- 
tually non-existent in this color range. 



4 SPATIAL DISTRIBUTION 

In Figure [S] we show the spatial distribution of each individ- 
ual stellar populations - selected using the boxes shown in 
the left panels of Figure [4]- and the corresponding contours. 
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Figure 6. Zoomed-in spatial distribution of BSG, RSG, and AGB stars in the Subaru FOV, and RGB from HST/ACS data. Several 
small holes in the stellar distributions are due to highly saturated stars and their associated bleed spikes. The inner Hi contour of 
Figure[5] as well as the outline of the combined footprint of the two ACS pointings, are shown. The cross and plus symbol in the bottom 
left panels represent the dynamical centre and the centre of the H I contour, respectively, to highlight the offset between the two. 



From top to bottom, the panels show the distribution of the 
MS+BSG, RSG, AGB, and RGB. To serve as gui des, in 
each panel we also plotted two Hi contours from the BC02 
study, representing densities of 0.1 and 2xl0 20 atoms cm -2 : 
the first one is their lowest density contour, while the second 
is the innermost contour that includes most of the observed 
MS+BSG stars belonging to HoII. We find that the shape of 
the latter contour is close to an ellipse with radius R=6.0', 
b/a=0.90, and position angle of 27°. These are the contours 
we used to separate HoII and field stars (see Section [3]). 

The top panels show that blue stars younger than about 
160 Myr old are confined to regions where the Hi density 
is higher than about 2xl0 20 atoms cm~ 2 . The very few ob- 
jects outside the contours are distributed uniformly over the 
FOV, suggesting that they are either unresolved galaxies 
with similar colors, or foreground blue HB, blue staggler, 
or white dwarf stars. The outermost MS+BSG contour fol- 
lows very closely the inner H I contour, as expected if mas- 
sive star formation only occurs above a certain gas den- 
sity threshold. The other panels show that contamination 
by field stars and unresolved galaxies increases significantly 
at redder color and fainter magnitudes. Nevertheless, these 
maps do not show significant stellar concentrations outside 
of the central contour, which suggests that most or all of 
the sources outside this area are either foreground stars or 
unresolved galaxies. 

Interestingly, in Figure [5] the concentration seems to 
increase from BSG to RSG stars, and from RSG to AGB 
stars, the populations being more compact with increas- 
ing age. The RGB stars, close to our detection limit, are 



the most affected by crowding which explains the hole in 
the distribution at the centre of HoII. In addition, some of 
the stars nagged as RGB may actually be RSG or AGB 
stars that were shifted to the RGB box due to photomet- 
ric errors. Therefore, since the RGB distribution is not clear 
from our Subaru data, we checked the distribution using the 
HST/ACS photometry of HoII described in Section |OJ In 
the bottom left panel of Figure [15] we show the distribution 
of RGB stars in the top 1.5 mag of the RGB from the HST 
dataset. The MS+BSG, RSG, and AGB maps from the Sub- 
aru data are shown in the other panels. In each panel, the 
black line outlines the combined footprint of the two ACS 
pointings to facilitate comparison. 

The difference between the distributions of the 
MS+BSG and RGB samples is striking: despite the small 
area covered by the HST data, it is obvious that the old, 
RGB stars are significantly more c oncentrated than th e 
younger stars, as already noted by iBastian et al.l |201l| ). 
The former are also distributed in a roughly circular dis- 
tribution, whereas the latter present a very irregular distri- 
bution. Within the ACS footprint, the density of RGB stars 
almost vanishes close to the edges while the density of blue 
stars is roughly uniform. We note that the distributions of 
Subaru and HST/ACS blue stars are very similar in the area 
in common. Using only the ACS data to avoid the uncertain- 
ties due to incompleteness and different spatial coverage, we 
find exponential scalelengths for the MS+BSG, RSG, AGB 
and RGB samples of 3.1' ±0.5', 1.7' ±0.3', 0.80' ±0.07', and 
0.76' ±0.04', respectively. Given the small FOV, we did not 
apply a background correction to these profiles. 
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Figure 7. Left: spatial distribution of MS+BSG stars overlaid on the Hi contours of iBCOsl Right: Surface density of MS+BSG stars as 
a func tion of H I densit y in square bins of 37.5" on a side. The solid line represent the typical Kennicutt-Schmidt relation with a slope 
of 1.4 l lKennicutdll998h . 



Figure [6] also shows that the blue stars, besides having 
a larger spatial extent, also seem to delineate short spiral 
arms. This morphology is not due to highly saturated stars 
or artifacts due to image defects, and completeness down 
to Io =24.5 is ~100 per cent at large radii. These arms 
follow the sa me counter-clockwise orientation as the H I arms 
observed bv lBureau et al.l (|2004l ) in the northwest of HoII, 
which further suggests they are real. The lack of similar 
arms in the distribution of AGB stars suggests these are 
likely transient features. 

The good agreement between the distribution of young 
stars and the H I gas is further illustrated in Figure [7] The 
left panel presents the dist ributio n of MS+BSG stars over- 
laid on the H I countours of lBC02l . and shows that these are 
usually located in the higher column density areas. In the 
right panel, we plot the number of MS+BSG stars in bins of 
37.5" on a side as a function of the H I density in the same 
bins. Note th at for this pan e l we h ave used the higher reso- 
lution map of IWalter et al.l l|2008l . natural weighting, beam 
size ~ 13") to better sample the scales of our stellar distri- 
bution. It shows that the distributions of young stars and 
Hi gas are correlated. A few points at high gas density and 
low star counts are artifacts due to the holes left by very 
saturated stars and their bleed spikes. On the other hand, 
the datapoints with high stellar and low gas density are real: 
most of them a re located in t he m ain H I hole at £ ~ 0' and 
r\ ~ —2', where IWeisz et all (|2009h found elevated star for- 
mation about 50 Myr ago, but only low levels in the past 
10-20 Myr. Interestingly, even though we use the number 
of MS+BSG stars (i.e. younger than about 150 Myr) as a 
rough proxy for the star formation rate, we find that the line 
representing the typical Ken nicutt-Schmidt relation with a 
slope of 1.4 (lKennicuttlll998l ) provides a good match to our 
data. 



Another interesting characteristic in Figure [6] is the 
apparent offset between the centres of the MS+BSG and 
HST/RGB star distributions, shown in the bottom left panel 
as the open circle and cross, respectively. We note that we 
applied a small shift (<4") to the astrometry of the HST 
data to correct for a slight mismatch between the two cat- 
alogues, so the offset is not due to inaccura te astrometry. 
We fin d that the dynamical centre listed in IStewart et al.l 
l|200d . 0:2000 = 8 h 19 m 05.6 s , <5 2 ooo = 70°43'25") provides an 
excellent fit to the distribution of RGB stars. On the other 
hand, the centre of the inner Hi contour, which follows very 
closely the distribution of MS+BSG stars, is offset from the 
dynamical centre by 30" (0.49 kpc) east and 7" (0.1 kpc) 
north. While the origin of this shift is unclear, it is reason- 
able to suspect that it is related to the processes that shaped 
the outer H I envelope. 



5 RADIAL PROFILE 

In order to obtain the radial density profile of HoII, one 
needs a reliable estimate of its ellipticity, position angle, 
and location of its centre. We use the dynamical cen- 
tre described in the previous Section, which is also the 
centre of the intermediate-age and old stellar populations 
(>1.5 Gyr old). The large uncertainty on the inclina- 
tion of HoII (IBC02I: Ide Blok et al 20081 : IWalter et al. 20081 : 



lOh et al. 201 it iGentile et al. 20121 )" prevents us from prop- 
erly constraining the morphology of the galaxy. However, 
given the very low ellipticity of the inner Hi contours 
(1— b/a~0.1) and circular distribution of the RGB stars, we 
can assume circular symmetry for the purpose of calculating 
the radial profiles. 

While resolved star counts are the optimal way to probe 
the very low surface brightness structure of galaxies (e.g. 
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iBarker et a,l.|[2009h . they are not ideal at small radii due to 
crowding and incompleteness. As described in Section 12.31 
crowding is particularly severe in the central part of HoII. 
However, it is possible to use diffuse light in this area and 
combine the contribution of both to obtain a composite pro- 
file extending over the whole radius. 

For the diffuse light surface-brightness profile, we used 
the median pixel value in concentric circles, after masking 
the saturated stars and their associated bleed spikes, as well 
as bright background galaxies. The sky value was estimated 
from the mode of the pixel value distribution of the en- 
tire image. The pixel values were then converted to mag- 
nitudes per arcsec 2 , calibrated using the same zero-point as 
for the stellar ph otometry, and corrected for reddenin g using 
E(B-V)=0.032 i|Schlegel. Finkbeiner. fc Davisll 19981 ). 

The same concentric circles were used for the star count 
profile, in which we counted the number of stars brighter 
than 7o = 24.5. This limit was chosen as the best com- 
promise between sufficient number statistics and reasonable 
completeness correction. The background contaminant level 
was estimated from the density of point sources in a wide 
circular annulus in a region where HoII's radial profile is 
flat within the uncertainties (7.5'^ r ^ f0.5'). For the area 
covered by the HST data, we applied the correction for com- 
pleteness determined in Section |2. 31 As for the diffuse light, 
the radial density profile E(r) was then converted to a mag- 
nitude scale using the relation fi(r) = —2.5 log E(r) + ZP. 
However, here the zero-point ZP was estimated by matching 
the overlapping region of the stellar and diffuse light profiles 
(2.5'< r < 4'). 

The resulting profiles are presented in Figure [8] the dif- 
fuse light and stellar surface-brightness profiles are shown as 
the gray filled circles and black filled stars, respectively. The 
open stars show the completeness-corrected star count pro- 
file. Although the completeness-corrected star counts and 
the diffuse light give reassuringly similar information in the 
inner regions, the star count data allow the profile to be 
extended to R ~ 7' where fj,y ~ 32 mag arcsec -2 . 

The surface brightness is well-described by an inner core 
(to ~1 .5') and a smooth outer decline. We find that an EFF 
l Elson. Fall, fc Freemanllf987h profile provides the best de- 
scription over the entire radial range, while an exponential 
fit also works b eyon d the inne r core . On the other hand, 
the iKind (I1962T ) and iPlummed (| 19 111 ) profiles (not shown) 
tend to under- and overestimate the stellar density beyond 
~3', respectively. The exponential fit to the profile between 
R=1.5 and 7' is shown as the dashed line, and has a scale- 
length of 0.70' ± 0.01', or 0.69 ± 0.01 kpc at the distance of 
HoII. The EFF fit (solid line) leads to a half-light radius of 
1.41' ± 0.04' (i.e. 1.39 ± 0.04 kpc) and absolute magnitude 
My = -16.3. 

Our scal elength is slightly smaller than the estimates 
from lSwatersI (/is = 0.99' a nd Kr = 0.86' after correcting for 
the difference in ellipticitv; 1 1999h . The discrepancy between 
the two estimates is most likely due to a combination of the 
shallower depth of their data ([ir < 27.5) and the steep- 
ening of our profile at larger radii. On the oth er hand, the 
radial profile of HoII shown in lOh et all (|201ll ) extends fur- 
ther out than the one presented here (out to R~9 kpc). We 
believe their extent may be overestimated as a consequence 
of using a relatively large ellipticity (b/a ~ 0.66), and a po- 
sition angle determined from the Hi data (PA=175°) which 
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Figure 8. Background-subtracted surface brightness profiles 
from the diffuse light (gray filled circles), resolved star counts 
(filled stars), and completeness corrected star counts ( open stars) . 
The open circles represent the H I density profile, from lBCOa . The 
dashed and solid lines are exponential and EFF profiles fit to the 
composite surface brightness profile (see text). 

is offset by abo ut 35° from th e position angle of the optical 
data (PA=30°; ISwaterdll999h . 

We thus find that the stellar component of HoII is sig- 
nificantly more compact than the H I cloud in which it is 
embedded, which extends out to ~ 16' (BC02). Even if some 
stars belonging to HoII are present beyond R ~ 7', their con- 
tribution to the galaxy luminosity is likely to be negligible. 
Assuming an exponential profile, such a population would 
represent less than 0.1 per cent of the total galaxy light. 



6 DISCUSSION 
6.1 Global Structure 

Using the resolved stars counts in the outer regions of HoII 
revealed an extended structure that is an order of magnitude 
fainter than what can be reached with diffuse light. Outside 
the central 2' area, the profile is well-fit by an exponen- 
tial model out to R=7' with a scalelength of 0.69 ±0.01 kpc. 
However, we have also found that the scalelength is different 
for each stellar population, in the sense that it decreases with 
increasing age. This is contrary to what is usu ally observed 
in dw arf galaxies in the nearby Universe (e.g. IZhang et al.l 
2012), where the region of active star formation has been 
shrinking with time. This implies that HoII would have been 
a more compact galaxy in the past. The very low density 
of RGB stars in the outskirts of the HST fields compared 
to that of the younger stars can be interpreted as a signifi- 
cantly lower SFR in these parts >1.5 Gyr ago (the minimum 
age of stars in the top magnitude of the RGB) than in the 
past ~160 Myr. This is in good agreement with the SFHs 
obtained from these HST d ata, which show a hig her SFR 
starting about 300 M yr ago (|Dalcanton et al]|2012 . see also 
iMcQuinn et al.ll2010l ). as well as a "drama tic rise in the SFR 
over the past 50 Myr" (|Weisz et alj|2008h . 
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Figure 9. Distribution of gal axies in the M81 group in standard coordinates, where the symbol size is proportional to the apparent 
integrated B-band magnitude (Karachent sev et alj |2002l. HoII is shown as the black circle at (0,0), and the main galaxies of the group 
are labeled. The outermost H I contour of HoII from BC02 is also shown, scaled up in size by a factor of 2. 



Given the large amount and relativ ely hig h density of 
the neutral hydrogen surrounding HoII (|BC02T ). it is some- 
what surprising that star formation was not more vigorous 
prior to about 300 Myr. This could suggest that the ex- 
tended Hi has recently been acquired or that it previously 
had too low a density for star formation to be significant. 
The process that shaped the spectacular H I envelope could 
have also increased the gas pressure in the inner regions of 
the galaxy and triggered the large scale star formation that 
we observe today. Two potential mechanisms to explain this 
are ram pressure by a hot IGM and a tidal interaction with 
a companion galaxy, which we discuss and compare below. 



6.2 Origin of the Hi morphology? 

HoII is embedded in a massive H I cloud that is significantly 
more extended than the optical counterpart. The remark- 
able morphology of the cloud in the low density outskirts, 
compressed on one side with a cometary appearance on the 
opposite side (see Figure [T}, is very suggestive of ram pres- 
sure stripping due to the presence of a hot IGM. While 
there are no X-ray observations to date confirming the pres- 
ence of diffuse hot gas in the vicinity of HoII, such observa- 
tions are available for other m uch smaller, 'poor' groups (e.g. 
IZabludoff fc Mulchaevlll998l ) in which the hot IGM density 
can be as high as ~ 5 x 10 -4 atoms cm -3 l|Mulchaev et al.l 
Il993h ISunll201^ . i.e. about two orders of magnitude higher 
than the minimum IGM density necessary to strip the inter- 
stell ar medium (ISM) of a low mass dwarf galaxy like HoII 
(e.g. |BC02| ; iMcConnachie et al. 20071 ). 

However, the disturbed appearance of the gas cloud may 
also be the result of gravitational interactions. While HoII 
is relatively isolated from the massive galaxies of the M81 
group (see Figure [9}, it has at least two companion dwarf 
galaxies within a projected distance of ~100 kpc which may 
have tidally affected the neutral hydrogen cloud. BC02 dis- 
cussed the possible origins of the Hi morphology and con- 



cluded that the ram pressure and tidal interaction scenarios 
were both plausible given the available information. Their 
follow-up study, despite the significant improvement in both 
coverage and sensitivi ty, was not sufficien t to rule-out one 
or the other scenario (|Bureau et al.ll2004l ). They suggested 
that deep stellar photometry could resolve the issue by re- 
vealing a stellar counterpart to the northwest H I structure, 
since ram pressure is not expected to have an effect on stars. 
Here we summarize the current arguments in favor of each 
scenario and discuss whether our deep photometry may help 
elucidate which is the main mechanism at play. 

The deeper Hi observations of lBureau et al.l (|2004l ) re- 
vealed that the low surface brightness component in the 
northwest of HoII can be resolved in two or three arms, 
which are often ass ociated with tidal interaction events (e.g. 
iDobbs et al1l2010l ). As shown in Section H the MS+BSG 
star distribution appears slightly distorted too, with spi- 
ral arm-like features south and west of HoII. The stel- 
lar arms also follow the same counter-clockw i se or ienta- 
tion as the Hi arms. However, iBureau et all l|2004t ) also 
noted that the arm kinematics followed the regular rota- 
tion of the inner gas, contrary to what would be expected if 
they were torn from the main body by tidal interactions. 
In addition, while the young stars do seem to trace spi- 
ral arms, the older populations have a much more regu- 
lar circular distribution. Tidal forces would affect all pop- 
ulations equally and the resulting stellar arms and/or tails 
would be made of stars of all ages. This suggests that most 
of the recent SF occured in the denser Hi arms, result- 
ing in a similar arm-like distribution of the young stars, 
rather than stellar tails due to tidal forces. The Hi arms, 
in turn, could be the consequence of ram pressure, as mod- 
els which include gas cooling have shown that it can pro- 
duce spiral arms (e.g. ISchulz fc Struckll200ll ; IVollmerl l 2003l ; 
iMapelli, Moore, fc Bland-Hawthornll200Sl ). 

Further in s ight c ome s from the detail ed SFHs of HoII 
bv lWeisz et~ai] i|200Sl ) and lDalcanton et al.l (|2012T ) , which re- 
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veal a significant enhancement in the last few hundred mil- 
lion years. Could this enhancement also be a consequence 
of the process that led to the striking morphology of the 
H I cloud? Strong bursts of star formation are often caused 
by an external trigger such as an accretion event or in- 
teraction with a nearby galaxy (e.g. iKennicutt et al.l 1 19871 ; 
Bernard et ai] |2012l ; ICignoni et alj|2012h . Indeed, as shown 
in Figure [9] HoII has two close companion dwarf galax- 
ies: M81dwA (=Kar52) located at a projected distance of 
~30 kpc to the North East, and UGC4483, ~110 kpc to the 
South East). The former, in particular, is sufficiently close 
that the timescale of a past interact ion is roughly compatibl e 
with the enhanced SFR found by iDalcanton et al.l |2012t ). 
Assuming a relative velocity of 100 km s -1 (the differenc e 
of radial velocities is ~45 km s~ 1 : lKarachentsev et al.l [2002). 
means that M81dwA could have passed close to HoII about 
300 Myr ago, in good agreement with the beginning of the 
SFR enhancement. Both companion dwarfs also appear to 
have been experiencing higher star formatio n in the past few 
hundred million years tha n at older epochs (|McQuinn et al.l 
l2010l : IWarren et alJlioTlh . 

However, tidal interaction is not the only process 
that can lead to enhanced star formation. In fact, hy- 
drodynamical and N-body simulations of the interaction 
between the hot IGM and the ISM of a galaxy have 
shown that ram pressure can also result in increased 
star formation in the inner regions. The Shockwaves re- 
sulting from the gas collision can increase significantly 
the central gas surface density and lead to the collapse 
of m olecular clouds, thus enhanc i ng the star formation 
rate dQuilis, Moore fc Bowerl |2000|; ISchulz fc Struck) l200ll ; 



IVollmer et al.ll200ll ; lMarcolini. Brighenti, fc D'Ercolelf2003f ). 

The fact that the spatial distribution of young stars closely 
follows that of the H I contours strongly supports this inter- 
pretation. Such an enhancement of the star formation as a 
consequence of the interaction between ISM and IGM has 
also been observed in various galaxies as asymmetric Ha en- 
hancements, usu ally located along the leadin g edge of the 
galax y disc (e.g. iKoopmann fc Kennevl |2004| ; ICrowl et al.l 
l2005l ). 

Therefore, while the morphology of the Hi cloud and 
the observed enhancement of the recent star formation rate 
cannot help discriminate between the phenomena at play, 
the regular circular distribution of the intermediate-age and 
old stars, the regular rotation of the Hi arms, the undis- 
turbed H 1 distribution of the two nearest dwarf compan- 
ions a nd the lack of Hi bridges/filaments between them and 
HoII l|Bureau et al.ll2004l ; IChvnoweth et al.ll2009l ; lOtt et all 
2012), all support the ram pressure stripping scenario. In 
addition, we do not find a stellar counterpart to the Hi 
cloud beyond R~7' where it becomes distorted and forms 
arms/tails, whereas tidal forces would affect gas and stars 
equally. Therefore, our data strongly suggest that ram pres- 
sure stripping is the main process responsible for the swept- 
back appearance of the Hi cloud. 



extent of the galaxy, including the vast H I cloud, and allow 
us to perform photometry of individual stars down to / ~ 
25.2, i.e. about 1.5 mag below the tip of the RGB. 

The deep CMDs reveal the presence of stellar popu- 
lations of all ages, from a few Myr old (MS+BSG, RSG) 
to several Gyr old (RGB). While in most dwarf galax- 
ies in the Local Universe the younger stars are found to 
be more centrally conc entrated than the older populations 
(e.g. IZhang et al]|2012j '). we find that in HoII the old RGB 
stars are significantly more concentrated than the young 
MS+BSG stars. Indeed, we find that the exponential scale- 
length for the young MS+BSG population is much larger 
than that of the RGB component (2.8'±0.4' vs. 0.76'±0.04', 
respectively). We speculate that the Shockwave due to ram 
pressure increased the gas density in the central part of the 
Hi cloud and triggered star formation on large scales. 

Our Subaru data enable us to construct a composite 
surface brightness profile for HoII by combining diffuse light 
in the central region with star counts at large radii. This 
profile is one of the deepest yet published for any galaxy, 
extending from the centre out to R ~ 7' where /j,v = 32 mag 
arcsec" 2 . Fitting an exponential profile to the outer regions 
gives a (projected) scalelength of 0.70' ±0.01', corresponding 
to 0.69 ± 0.01 kpc at the distance of HoII. 

Finally, we discuss the properties of the resolved stel- 
lar populations in the context of the morphology of the large 
H 1 cloud in order to understand the origin of its swept-back, 
cometary appearance. Previous studies based on diffuse op- 
tical light or 21 cm data could not definitively determine 
whether the cloud shape was due to ram pressure from a 
hot IGM or to a tidal interaction with a nearby companion 
galaxy. Our deep photometry shows that the intermediate- 
age and old stars have a regular circular distribution and 
show no sign of tidal tails/streams. In addition, we find that 
there are very few, if any, HoII stars beyond R~ 7' where 
the H 1 becomes distorted. Since tidal forces would affect gas 
and stars equally, our data strongly suggest that the spec- 
tacular morphology of the H I cloud is due to ram pressure. 
The detection of significant amount of diffuse hot gas in the 
vicinity of HoII would further verify this. 



ACKNOWLEDGMENTS 

We are very grateful to the anonymous referee for a prompt 
report that helped us improve the manuscript, and would 
like to thank M. Bureau for providing the Hi density map 
of HoII. Support for this work was provided by a rolling 
grant from the Science and Technology Facilities Coun- 
cil. We acknowledge the usage of the HyperLeda database 
(http://leda.univ-lyonl.fr). This research has made use of 
THIN GS, 'The Hi Nearby Galaxy Survey' (| Walter et all 
|200S| ). and the NASA/IPAC Infrared Science Archive, which 
is operated by the Jet Propulsion Laboratory, California 
Institute of Technology, under contract with the National 
Aeronautics and Space Administration. 



7 SUMMARY AND CONCLUSIONS 

We have carried out a wide-field survey of the M81 group 
dwarf galaxy Holmberg II based on deep Subaru/Suprime- 
Cam imaging in V and I. These observations cover the whole 
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